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The observation of double parton collisions by CDF has provided the first 
^^ • direct information on the structure of the proton in transverse space. The 

[T • actual quantity which has been measured is the 'effective cross section' (Jeff: 

vQ ■ which is related to the transverse size of the region where hard interactions 

CN . are localized. The actual value which has been measured is sizably smaller 

than naively expected and it is an indication of important correlation effects 
^^ ■ in the many-body parton distribution of the proton. We discuss the problem 

t~^ ■ pointing out a possible source of correlations in the proton structure, which 

^ . could have a significant effect on the value of cJe// . 
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I. INTRODUCTION 

While parton distributions represent the relevant information on the non perturbative 
input for most of the large pt processes at present energy, parton distributions do not exhaust 
the information on the hadron structure. In fact the information in the parton distributions 
corresponds to the average number of partons of a given kind and with a given momentum 
fraction x which is seen when probing the hadron with the resolution Q^. In this respect 
a qualitative change has occurred when it has been shown that the observation of hadronic 
interactions with multiple parton collisions is experimentally feasible. In a multiparton 
collision event different pairs of partons interact independently at different points in the 
transverse plane [|^ and the process, as a consequence, depends in a direct way on the actual 
distribution of the interacting partonic matter in transverse space. The simplest case is the 
double parton scattering. The non perturbative input to the process is the double parton 
distribution D2{x,x';h) depending on the momentum fractions of the interacting partons 
and on their relative distance in the transverse plane b. If all partons are uncorrelated and 
the dependence on the different degrees of freedom is factorized, one may write D2{x, x'; b) = 
geff{x)geff{x')F(h) where geff{x) is the usual effective parton distribution, namely the gluon 
plus 4/9 of the quark and anti-quark parton distributions. When the interacting hadrons 
are identical F{h) is equal to the overlap in transverse plane of the matter distribution of 
the two hadrons (normalized to one) as a function of b, which now represents the impact 
parameter. The effective cross section aejf is introduced by the inclusive double scattering 
process, which is proportional to l/cTeff = J (PbF'^{h) |jl|]. The effective cross section is then 
a well defined property of the hadronic interaction and, at least if the simplest possibility 
for the hadronic structure is realized, it is both energy and cutoff independent. Analogously 
the dimensional scale factors, which one may introduce in relation to triple, quadruple etc. 
partonic collisions, are energy and cutoff independent properties of the interaction. 

Initially the search of double parton collisions have been sparse and not very consistent 
0. CDF however has recently claimed the observation of a large number of double parton 
scatterings [§]. The measured value of the effective cross section, a^fj = 14.5 ± l.Tlll'Imb, is 
sizably smaller than expected naively and it is an indication of important correlation effects 
in the hadron structure [Q . Correlations on the other hand are a manifestation of the links 
between the different constituents of the hadronic bound state and, in fact, it is precisely 
this sort of connections which one would like to learn as a result of experiments probing the 
hadron structure. 

In the present paper we point out a rather natural possible source of correlations, which is 
able to give rise to a sizable reduction of the value of the effective cross section, as compared 
to the most naive expectation. In fact we show that, when linking the transverse size of the 
gluon and sea distributions to the actual configuration of the valence, the expected value of 
(Te// is sizably decreased. In the next paragraph we recall the main features of the multiple 
partonic collision processes in the simplest uncorrelated case, in the following section we 
correlate the distributions of gluons and sea quarks with the configuration of the valence 



and, in the final paragraph, we discuss our results. 

II. MULTIPLE PARTON COLLISIONS 

The non-perturbative input to the multiple parton collision processes is the many body 
parton distribution. In most cases the simplest uncorrelated case, with factorized dependence 
on the actual degrees of freedom, is considered and the probability distribution is expressed 
as a poissonian. The probability of a configuration with n partons with fractional momenta 
xi, . . . , Xn and with transverse coordinates bi, . . . , b„ is then expressed as 

r(xi, bi, . . .x„, b„) = ''- j — — ^exp { — / -D(x, h)dxd%\ (1) 

where -D(x, b) = geff{x)f{h) is the average number of partons with momentum fraction x 
and with transverse coordinate b. The function geff{x) is the effective parton distribution 
while /(b) represents the partonic matter distribution in transverse space. The many body 
parton distribution is an infrared divergent quantity, one needs therefore to introduce an in- 
frared cut-off. A natural cutoff in this context is the lower value of momentum transfer which 
allows a scattered parton to be recognized as jet in the final state. Given the many-body 
parton distribution one can work out all possible multi-parton interactions. If rescatterings 
are neglected, namely if every parton is allowed to interact with large momentum transfer 
(namely larger than the infrared cutoff) at most once, one may write a simplest analytic 
expression for the hard cross section an, corresponding to the cross section counting all 
inelastic hadronic events with al least one hard partonic interaction 0: 
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here (3 is the impact parameter of the hadronic collision, F{[3) = J (fbf(h)f(h — (3) and as 
is the usual expression for the integrated inclusive cross section to produce jets, namely the 
convolution of the parton distributions and of the partonic cross section. The unitarized 
expression of an in Eq.^, which takes into account the possibility of many parton-parton 
interactions in an overall hadron-hadron collision, is well behaved in the infrared region and 
it corresponds to a poissonian distribution of multiple parton collisions at a fixed value of 
the impact parameter of the hadronic interaction. If one works out in Eq.^ the average 
number of partonic collisions one obtains 

{n)aH = Jd'(3asF{(3) = as (3) 

so that as represents the inclusive cross section normalized with the multiplicity of partonic 
interactions. The inclusive cross section for a double parton scattering ao, normalized in an 
analogous way with the multiplicity of double parton collisions, {n{n — l)/2)aH, is given by 

^'"^'"~^^^ ^H = lJd^P4F\p) = an (4) 



One may then introduce the effective cross section: 

1 2 

which is therefore expressed in terms of the overlap of the partonic matter distribution of 
the two interacting hadrons as 

In fig.|I|two different analytic forms for /(b) are compared: a gaussian and the projection 
of a sphere in transverse space. The radius of the distributions has been fixed in such a way 
that in both cases the RMS radius has the same value that, following CDF 0, has been set 
equal to 0.56fm. For the cutoff we have used 5GeV and, in order to be able to reproduce the 
value of the integrated hard cross section for producing minijets measured by UAl |^ , we 
have multiplied the partonic cross section, computed in the perturbative QCD parton model, 
by the appropriate fc-factor. Both choices give a similar qualitative result: the effective cross 
section is constant with energy (and cutoff independent) while the 'total' hard cross section 
an grows rapidly as a function of the cm. energy. The values obtained for aeff are however 
too large, by roughly a factor two, as compared to the value quoted by experiment. One 
has to add that in the experimental analysis of CDF all events with triple parton collisions 
have been removed from the sample of inelastic events with double parton scatterings. The 
resulting value for the effective cross section, o"e//|cDF, is therefore larger with respect to the 
quantity discussed here above and usually considered in the literature. The disagreement 
with the most naive picture is therefore even stronger then apparent when comparing the 
result of the uncorrelated calculation with the quoted value of o"e//. 

In the usually considered picture of multiparton interactions just recalled all correlations 
are neglected and one may therefore claim that the experimental evidence is precisely that 
correlations play an important role in the many-body parton distribution of the hadron. In 
the next paragraph we propose therefore a slight modification to the picture, linking the 
transverse size of the gluon and sea distributions to the configuration in transverse space of 
valence quarks. 

III. A SIMPLEST MODEL FOR THE PARTONIC MATTER DISTRIBUTION IN 

THE PROTON 

Charged matter is distributed in the proton according to the charge form factor, which 
is well represented by an exponential expression in coordinate space. The information refers 
to a large extent to the distribution in space of the valence quarks, which can then be found 
in various different configurations in transverse space with a given probability distribution. 
Less is known on the distribution in space of the neutral matter component of the pro- 
ton structure, namely the gluons. In this respect one may consider two different extreme 
possibilities: 



• the distribution in space of sea quarks and gluons has no relation with the distribution 
of valence quarks, 

• or, rather, its transverse size is linked closely to the actual configuration in space of 
the valence quarks. 

The no correlation hypothesis, which has been ruled out by the measure of <y(.ff, would 
imply the first possibility. We try therefore to work out here a simplest model where the 
distribution of the whole partonic matter in the hadron is driven by the actual configuration 
of valence quarks. 

The parton distribution in Eq.^ can be modified by separating the valence quarks from 
gluons and sea quarks and by correlating the transverse size of the distributions of sea quarks 
and gluons with the transverse size of the actual configuration of the valence. If however 
one simply rescales the distribution function in transverse space /(b) by rescaling the radius 
while keeping the normalization constant, one is imposing a conservation constraint in the 
number of gluons. The number is in fact the same in each configuration of the proton, both 
when it is squeezed to a small transverse size and when it is expanded to a relatively large 
dimension. The picture therefore is not consistent with the common believe that the energy 
to the gluon field grows basically because of the growth of the distance between the valence 
quarks. It seems therefore more reasonable to remove the constraint on normalization and 
to keep rather fixed the maximum value of /(b) in all configurations. We modify therefore 
the simplest poissonian expression in Eq.|l| and we express the probability distribution T as 
follows: 



r(xi, bi, . . . Xn, b„) = (fiBo, B)q^{Xi)q^{X2)qv{X: 
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To simplify the notation we have not written explicitly the dependence of F on the coordi- 
nates of the valence quarks, X and B. The dependence of F on the momentum fraction of 
the valence quarks X is factorized and qv{X) is the usual distribution of valence quarks as 
a function of the momentum fraction X. The function f{b, hi) represents the distribution 
of gluons and sea quarks in transverse space. It is a function of the transverse coordinate of 
the considered parton, bj, and it depends on the scale factor b whose value is given by the 
actual configuration of the valence in transverse space. The transverse coordinates of the 
three valence quarks are 



Bi = 


-B, + B 


B2 = 


iB,-B 


B3 = 


-Bd 



(8) 

The dependence on the transverse coordinates of the valence is given by ip(Bi:),B), that 
is the integral on the longitudinal coordinates Zo, Z of 0(Rd, R), representing the valence 
structure of the proton in coordinate space. Explicitly we use the exponential form 



A^ A^ 

(RB,R) = T^exp{-(Az?i?D + Ai?)} (9) 



where 
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and J {r"^) = O.Slfm is the proton charge radius. 

In a given configuration of the valence, the average number of gluons and sea quarks is 
not equal to the overall average number g{x) (with g{x) we indicate here the sum of gluon 
and sea quark distributions). It is rather equal to g{x)-^; where (6^) is the average of 6^ 
with the probability distribution of the valence. To make a definite choice we take b = Bo 
in such a way that 

{I?) = {BD = I dRndIiBl(l){Rn, R) (11) 

The density of gluons and sea quarks in transverse space is then constant in the middle 
of the proton in all different configurations of the valence quarks and it is equal to the 
value assumed in the average configuration. The actual number of gluons and sea quarks is 
therefore a function of the configuration taken by the valence, compact configurations giving 
rise to small numbers while more extended configurations giving rise to larger numbers. 

Summing over all possible probability configurations of gluons and sea quarks one ob- 
tains, from Eq.^, the probability distribution of valence in transverse space (p(Bd, B). The 
average number of gluons and sea quarks, with given momentum fraction x, is 

JdBndBdh^{BD,B)g{x)f{Bn,h)^=9i^) (12) 

while the average number of partonic collisions, non involving valence {n{B£), B, B'^,, B', (3)), 
for a given configuration of the valence and for a given value of the impact parameter /3, is 
expressed as 

{niBn,B,B'^,B',P))= fdhf{Bn,h-P)^fiB'^,h)^-^[asig + qs,g + qs)] (13) 

The average number of interactions of valence quark with gluons and sea is written in an anal- 
ogous way. If one integrates the average number of partonic collisions (n(B£), B, B^, B', (3)) 
over all configurations of the valence quarks of the two hadrons, with the corresponding 
weights, and on the hadronic impact parameter /3, one obtains the single scattering inclu- 
sive cross section as'- 

I dBodBdB'^dB'dPipiBD, B)if{B'^, B') (n(B^, B, B^, B', p)) = asig + g„ g + g.) (14) 



An analogous expressions for interactions of valence quarks with gluons and sea quarks is 
readily written. 

All average quantities are therefore the well known ones of the perturbative QCD-parton 
model. The inclusion of the transverse degrees of freedom in the relations above allows one 
to write down promptly all expressions corresponding to the various multiparton collision 
processes. For the double parton scattering case one has: 

an = ^y'ciB^dBrfB'^dB'd/3(^(Bz5,B)(^(B^,B')[(n(B,5,B,B'^,B',/3))]' (15) 

and for the hard cross section 

an = f dBDdBdBodB'd(3ip{BD,B)ip{B'^,B')[l - exp[ - (^(Bd, B, B'^, B',;^))}] (16) 

The effective cross section is obtained, as in the uncorrelated case, from an by using Eq.||. 

IV. DISCUSSION 

The values of an and of CTeff, derived from the correlated parton distribution described in 
the previous paragraph, are plotted in fig.H, where two different analytic expressions for /(b) 
are considered. In figj^a /(b) is the projection of a sharp-edged sphere in transverse space, 
while in fig.^fe /(b) is a gaussian. In each figure we draw the inclusive cross section as, which 
is the fast growing short-dashed curve, the hard cross section an, continuous curves, and 
the effective cross section <Jeff, long-dashed curves. Both for an and for o"e// we draw two 
different curves. The higher curve refers to the case where the gluon and sea distributions 
are kept fixed and equal to the average configuration, irrespectively of the configuration of 
the valence, while the lower curve refers to the case where the configuration of gluons and sea 
are correlated in radius with the configuration of the valence, as described in the previous 
paragraph. The lower threshold for the transverse momentum of the produced jets has been 
put equal to 5GeV. In the correlated case, when considering a sphere for /(b), we choose 
for the radius of the sphere the value of Bd. The RMS of the radius of the sphere, averaged 



with the probability distribution of the valence, is then J {B'j^) = 0.66fm, namely it is equal 



to -i/2/3(r^), with J {r^) = O.Slfm, the RMS proton charge radius. In the gaussian case 
we fix the size of the distribution by requiring the same RMS value for the radius of the 
distribution as in the case of the sphere. As one may see in fig.0, while the value of o"e// is 
too large when the distribution of gluons and sea quarks is kept fixed, the value of aejf is 
sizably reduced when it is correlated with the distribution of the valence quarks. 

With our choice the value of the effective cross section turns out of the right size. Ob- 
viously the result could have been significantly different with a different, but still plausible 
choice of the (correlated) radius of the gluon and sea distributions. We do not claim there- 
fore that the simplest mechanism discussed here is the only solution to the problem posed 
by the smallness of the observed value of Ceff- Rather we point out that the source of corre- 
lations discussed here on one side is a minimal modification to the uncorrelated many-body 



parton distribution and, on the other, it looks as a natural possibility, which could explain 
a substantial amount of the difference between the observed value of o"e// and the result of 
the uncorrelated calculation. 

As shown in fig.^, our simplest model, while reducing sizeably the expectation for (Jeff, 
does not modify dramatically the expectation for an- It is therefore worth wile making 
predictions for other possible observables, in order to have an independent indication on the 
actual model. We have then estimated the triple parton scattering cross section. The triple 
scattering cross section, being proportional to o"|, introduces a new dimensional quantity 
other than the effective cross section Ceff- In the uncorrelated case the new dimensional 
quantity has a given value, proportional to a^j^j and the proportionality factor depends on 
the actual form of /(b). Also in the correlated case one may however write 

"T = J,^ (IT) 

The observation of the triple scattering parton process allows therefore one to measure the 
dimensionless quantity r, which, as (Te// is a well defined quantity, related to the geometrical 
properties of the interacting hadrons. The actual expectations for the r-factor according to 
the model discussed here are shown in the table. 



V. CONCLUSIONS 

The observation of multiple parton collisions allows the measure of a whole set of quan- 
tities which characterize the interaction and that are directly connected to the geometrical 
properties of the hadron structure. The first indication in this direction is the measure 
of the effective cross section of double parton collisions performed by CDF. The measured 
value of the effective cross section rules out the simplest uncorrelated picture of the many 
body parton distribution of the hadron. In this note we have shown that correlating the 
transverse dimension of the gluon and sea quark distributions to the transverse dimension 
of the valence one can obtain for the effective cross section a value much closer to the ex- 
perimental indication. In order to have the possibility to test the model, at least to some 
extent, we have then worked out our prediction for the factor r, characterizing the triple 
parton scattering cross section. 

The simplest expectation for (Xe// and, more in general, for the scale factors of the 
different multiple parton collision processes, is that they are cutoff and energy independent. 
The situation however is changed when more elaborate structures are considered. Also in 
the simplest model discussed here in fact aeff shows a slight energy and cutoff dependence. 
The origin is the following: In our case partons are organized in two different structures in 
transverse space, on one side one has valence quarks and on the other gluons and sea quarks. 
When changing the energy or the lower cutoff in pt one is varying the relative number of 
interacting sea quarks and gluons with respect to the valence quarks. This variation is then 



reflected in the relative weight of the rate of collisions of partons which have a different 
distributions in transverse space and, as a consequence, a^ff is slightly modified. 
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TABLES 





TABLE I. Factor r for the triple parton scattering process 


Sphere (no corr.) 


Gauss (no corr.) Sphere (with corr.) Gauss (with corr.) 


0.78 


0.74 0.5 0.46 
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FIGURES 
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FIG. 1. (J5, an and dg// in the simplest uncorrelated case. The two curves for an and o"e// 
refer to the two different choices made for /(b), a gaussian (higher curve) and a sphere (lower 
curve) . 
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FIG. 2. (T5, (T/f and (Te// when taking for /(r) a sphere a) or a gaussian 6). The lower curves 
for an and a^fj refer to the correlated many-body parton distribution, the higher curves to the 
uncorrelated ones. 
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